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The gating charge should not be estimated by fitting a two-state 
model to a Q-V curve 
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The voltage dependence of charges in voltage-sensitive proteins, typically displayed as charge versus voltage ( Q-V) 
curves, is often quantified by fitting it to a simple two-state Boltzmann function. This procedure overlooks the fact 
that the fitted parameters, including the total charge, may be incorrect if the charge is moving in multiple steps. 
We present here the derivation of a general formulation for Q-Vcurves from multistate sequential models, includ- 
ing the case of infinite number of states. We demonstrate that the commonly used method to estimate the charge 
per molecule using a simple Boltzmann fit is not only inadequate, but in most cases, it underestimates the moving 
charge times the fraction of the field. 
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Many ion channels, transporters, enzymes, receptors, 
and pumps are voltage dependent. This voltage depen- 
dence is the result of voltage-induced translocation of 
intrinsic charges that, in some way, affects the confor- 
mation of the molecule. The movement of such charges 
is manifested as a current that can be recorded under 
voltage clamp. The best-known examples of these cur- 
rents are "gating" currents in voltage-gated channels 
and "sensing" currents in voltage-sensitive phosphatases. 
The time integral of the gating or sensing current as a 
function of voltage (V) is the displaced charge Q_(V) , 
normally called the Q-Vcurve. 

It is important to estimate how much is the total 
amount of net charge per molecule (Q, !( «) that relocates 
within the electric field because it determines whether a 
small or a large change in voltage is necessary to affect 
the function of the protein. Most importantly, knowing 
Qjnax is critical if one wishes to correlate charge move- 
ment with structural changes in the protein. The charge 
is the time integral of the current, and it corresponds to 
the product of the actual moving charge times the frac- 
tion of the field it traverses. Therefore, correlating charge 
movement with structure requires knowledge of where 
the charged groups are located and the electric field 
profile. In recent papers by Chowdhury and Chanda 
(2012) and Sigg (2013), it was demonstrated that the 
total energy of activating the voltage sensor is equal to 
Qjnax Vm, where Vm is the median voltage of charge trans- 
fer, a value that is only equal to the half-point of activa- 
tion Vi/2 for symmetrical Q-Vcurves. Vm is easily estimated 
from the Q- Vcurve, but Qj nax must be obtained with other 
methods because, as we will show here, it is not directly 
derived from the Q-Vcurve in the general case. 

The typical methods used to estimate charge per 
molecule Q^ax include measurements of limiting slope 
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(Aimers, 1978) and the ratio of total charge divided by 
the number of molecules (Schoppa et al., 1992). The 
discussion on implementation, accuracy, and reliability 
of these methodologies has been addressed many 
times in the literature, and it will not be discussed here 
(see Sigg and Bezanilla, 1997). However, it is worth 
mentioning that these approaches tend to be techni- 
cally demanding, thus driving researchers to seek alter- 
native avenues toward estimating the total charge per 
molecule. Particularly, we will discuss here the use of a 
two-state Boltzmann distribution for this purpose. Our 
intention is to demonstrate that this commonly used 
method to estimate the charge per molecule is gener- 
ally incorrect and likely to give a lower bound of the 
moving charge times the fraction of the field. 

The two-state Boltzmann distribution describes a 
charged particle that can only be in one of two positions 
or states that we could call Si and S?- When the particle 
with charge Q, MI (in units of electronic charge) moves 
from Si to Sa, or vice versa, it does it in a single step. The 
average charge found in position S%, Q( V) , will depend on 
the energy difference between Si and &, and the charge of 
the particle. The equation that describes Q( V) is: 
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where V^is the potential at which the charge is equally 
distributed between Si and S2, and k and T are the 
Boltzmann constant and absolute temperature, respec- 
tively. The Q(V) is typically normalized by dividing Eq. 1 by 
the total charge Q MX . The resulting function is frequently 
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called a "single Boltzmann" in the literature and is used 
to fit normalized, experimentally obtained QV curves. 
The fit yields an apparent V1/2 ( V V2 ) and an apparent 
Qjviax ( CLax ) > an d this last value is then attributed to be 
the total charge moving (9 max . Indeed, this is correct but 
only for the case of a charge moving between two posi- 
tions in a single step. However, the value of Q_ max thus 
obtained does not represent the charge per molecule for 
the more general (and frequent) case when the charge 
moves in more than one step. 

To demonstrate the above statement and also esti- 
mate the possible error in using the fitted Q™* from Eq. 1, 
let us consider the case when the gating charge moves 
in a series of n steps between n + 1 states, each step with 
a fractional charge z, (in units of electronic charge eg) 
that will add up to the total charge Qj, la x- 
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The probability of being in each of the states Si is la- 
beled as P h and the equilibrium constant of each step is 
given by 



ft = exp 



z t (V-Vi) 
kT 



where z, is the charge (in units of e 0 ) of step i, and V, is the 
membrane potential that makes the equilibrium constant 
equal 1. In steady state, the solution of P, can be obtained 
by combining 



5±i 



■ p., i = 1... n and V P i = 1, 



giving 



m — 1 

i+zrift 



■ l..,n and P 1 



1 

n j 



i+irift 



We define the reaction coordinate along the moved 



charged ^as 



The Q-Vcurve is defined as 

n 
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Then, replacing P s yields 
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or written explicidy as a function of V: 
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Eq. 2 is a general solution of a sequential model with 
n + 1 states with arbitrary valences and Vis for each tran- 
sition. We can easily see that Eq. 2 has a very different 
form than Eq. 1 , except when there is only a single tran- 
sition (n= 1). In this latter case, Eq. 2 reduces to Eq. 1 
because Z; and V\ are equal to Q, MX and V1/2, respectively. 
For the more general situation where n > 1 , if one fits 
the Q(V) relation obeying Eq. 2 with Eq. 1, the fitted 
Q_ max value will not correspond to the sum of the z, val- 
ues (see examples below and Fig. 1). A simple way to 
visualize the discrepancy between the predicted value of 
Eqs. 1 and 2 is to compute the maximum slope of the 
Q- V curve. This can be done analytically assuming that 
Vi = V D for all transitions and that the total charge Qj nax is 
evenly divided among those transitions. The limit of the 
first derivative of the Q( V) with respect to devaluated at 
V= V Q is given by this equation: 



dV , 



\1nkT 



(3) 



From Eq. 3, it can be seen that the slope of the (9-Vcurve 
decreases with the number of transitions being maxi- 
mum and equal to Qjn ax /(4kT) when n = 1 (two states) 
and a minimum equal to Qj nra /(12kT) when n goes to in- 
finity, which is the continuous case (see next paragraph) . 

Infinite number of steps 

Eq. 2 can be generalized to the case where the charge 
moves continuously, corresponding to an infinite num- 
ber of steps. If we make 



-v 0 , 



then all p, t = p,, and we can write Eq. 2 as the normalized 
Q(V) in the limit when n goes to infinity: 
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Eq. 4 can also be written as 
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which is of the same form of the classical equation of 
paramagnetism (see Kittel, 2005). 

Examples 

We will illustrate now that data generated by Eq. 2 can 
be fitted quite well by Eq. 1, thus leading to an incorrect 
estimate of the total charge moved. Typically, the experi- 
mental value of the charge plotted is normalized to its 
maximum because there is no knowledge of the abso- 
lute amount of charge per molecule and the number of 
molecules. The normalized QV curve, Q lor , is obtained 
by dividing Q(V) by the sum of all the partial charges. 

Fig. 1 shows Qpor computed using Eq. 2 for one, two, 
three, four, and six transitions and for the continuous 
case using Eq. 5 (squares) with superimposed fits to a 
two-state Boltzmann distribution (Eq. 1, lines). The com- 
putations were done with equal charge in each step (for 
a total charge Qj„ ax = 4eo) and also the same K = — 25 mV 
value for all the steps. It is clear that fits are quite ac- 
ceptable for cases up to four transitions, but the fit sig- 
nificantly deviates in the continuous case. 

Considering that experimental data normally have 
significant scatter, it is then quite likely that the experi- 
menter will accept the single-transition fit even for cases 
where there are six or more transitions (see Fig. 1 ) . In 
general, the case up to four transitions will look as a 
very good fit, and the fitted Q mx value may be inaccu- 
rately taken and the total charge transported might be 



underestimated. To illustrate how bad the estimate can 
be for these cases, we have included as insets the fitted 
value of Qjmax for the cases presented in Fig. 1. It is clear 
that the estimated value can be as low as a fourth of the 
real total charge. The estimated value of V 1/2 is very 
close to the correct value for all cases, but we have only 
considered cases in which all V;'s are the same. 

It should be noted that if jm, of the rightmost transi- 
tion is heavily biased to the last state ( V, is very nega- 
tive) , then the Q max estimated by fitting a two-state model 
is much closer to the total gating charge. In a three-state 
model, it can be shown that the fitted value is exact 
when Vj— >°° and Vt— >— 00 because in that case, it con- 
verts into a two-state model. Although these values of V 
are unrealistic, the fitted value of Q m „ can be very close 
to the total charge when V 2 is much more negative than 
F; (that is, Vx » V 2 ) ■ On the other hand, If « V 2 , the 
Q- Vcurve will exhibit a plateau region and, as the differ- 
ence between Vj and V 2 decreases, the plateau becomes 
less obvious and the curve looks monotonic. These cases 
have been discussed in detail for the two-transition 
model in Lacroix et al. (2012). 

We conclude that it is not possible to estimate un- 
equivocally the gating charge per sensor from a "single- 
Boltzmann" fit to a Q-V curve of a charge moving in 
multiple transitions. The estimated Q mx value will be a 
low estimate of the gating charge Qp BX , except in the case 
of the two-state model or the case of a heavily biased late 
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Figure 1 . Examples of normalized QV curves for a Q lax = 4 computed with Eq. 2 for the cases of one, two, three, four, and six transitions 
and the continuous case using Eq. 5 (squares). All the QV curves were fitted with Eq. 1 (lines). The insets show the fitted valence ( Q max ) 
and half-point ( V 1/2 ) . 
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step, which are rare occurrences. It is then safer to call 
"apparent gating charge" the fitted Ci™* value of the 
single-Boltzmann fit. 

Addendum 

The most general case in which transitions between 
states include loops, branches, and steps can be derived 
directly from the partition function and follows the 
general thermodynamic treatment by Sigg and Bezanilla 
(1997), Chowdhury and Chanda (2012), and Sigg 
(2013). The reaction coordinate is the charge moving 
in the general case where it evolves from q= Oto q= 
by means of steps, loops, or branches. In that case, the 
partition function is given by 



z=Z ex P 
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kT 



(6) 



We can compute the mean gating charge, also called the 
Q-V curve, as 

Q{V ) = ( q ) = kT^ = kT^ = ^ } kT - (7) 
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The slope of the Q-Vis obtained by taking the derivative 
of <q> with respect to V. 
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Let us now consider the gating charge fluctuation. 
The charge fluctuation will depend on the number of 
possible conformations of the charge and is expected 
to be a maximum when there are only two possible 
charged states to dwell. As the number of intermediate 
states increases, the charge fluctuation decreases. Now, 
a measure of the charge fluctuation is given by the vari- 
ance of the gating charge, which can be computed from 
the partition function as: 



(A q *) = ( q >)-( q f= { kTf 
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But the variance (Eq. 9) is identical to the slope of Q(V) 
(Eq. 8) . This implies that the slope of the Q-Vis maximum 
when there are only two states. 
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